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Abstract -Two significant consequences of quantum fluctuations are entanglement and critical- 
ity. Entangled states may not be critical but a critical state shows signatures of universality in 
entanglement. A surprising result found here is that the entanglement entropy may become arbi- 
trarily large and negative near the dissociation of a bound pair of quantum particles. Although 
apparently counter-intuitive, it is shown to be consistent and essential for the phase transition, 
by mapping to a classical problem of DNA melting. We associate the entanglement entropy to a 
subextensive part of the entropy of DNA bubbles, which is responsible for melting. The absence 
of any extensivity requirement in time makes this negative entropy an inevitable consequence of 
quantum mechanics in continuum. Our results encompass quantum critical points and first-order 
transitions in general dimensions. 



C 

O 

CM 
> 

m 
(N 

d 



X 



Introduction. — Quantum entanglement [HI] is a 
fundamental feature of quantum mechanics, which says 
that performing a local measurement may instantaneously 
affect the outcome of local measurements far away. There 
is another feature of quantum mechanics where the zero- 
point fluctuations in the ground state may coherently add 
up to produce long-range correlations of local observablcs. 
This happens at quantum critical points (QCP), a point 
where the spectrum becomes gaplcss, obtained by tun- 
ing the parameters of the Hamiltonian. In both 
pure state cannot be written as the product of the wave- 
functions of the two distant parts, though states may be 
entangled without being critical. The ground-state energy 
may be non-analytic through a quantum phase transition 
(QPT) or through a quantum critical point. The wave 
function encodes not only this non-analyticity but also the 
special quantum correlations or quantum entanglement in- 
trinsic to the state. 

At or near a QCP, the signatures of its universality can 
therefore be found in the entanglement, a common mea- 
sure of which is the von Neumann entropy (5) [ME]- The 
exact results of this paper show that for a class of critical 
points, viz., the dissociation of a pair of particles in the 
unitarity limit of infinite scattering length, there is the 
possibility of a negatively diverging S. Although counter- 
intuitive, this is not an artifact. Analogous situation oc- 



curs in statistical mechanics for Gibbs entropy in canoni- 
cal ensemble for a gapless spectrum. As discussed below, 
the problem in hand involves a gapless entanglement spec- 
trum. The usual proof of the positivity of entanglement 
entropy is not applicable in case of continuous eigenvalues 
of the reduced density matrix. The negative entropy is es- 
sential for the criticality itself. Its importance is brought 
out via the mapping of the quantum problem to the equiv- 
alent classical statistical mechanical problem, the melting 
of a double-stranded DNA l9rfl2l. 



Entanglement entropy. — Recall the problem of a 
quantum particle of mass to in a three dimensional spher- 
ical potential well, 



V(r) = ~Vq for r < a, 
= for r > a, 



(1) 
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where r is the radial coordinate, a and Vq are the width 
and the depth of the potential well. What is special is that 
Vb > does not guarantee the existence of a bound state, 
unlike in one or two dimensions, or in classical mechanics. 
No bound state exists for u < u c where u = 2mVoa 2 /h 2 is 
the dimcnsionless parameter for the potential and u c cor- 
responds to a critical value of u. For simplicity, we take 
u ~ u c so that there is only one bound state. In this situa- 
tion energy \E\ itself is the gap in the spectrum. If we tune 
u to get a state with zero energy [E = 0), then at that 
energy in d = 3 the wave function tp(r) ~ l/r which is like 
a non-normalizable critical state. Like a bound state the 
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Fig. 1: (a) Gap A in the energy spectrum. The shaded region 
is the continuum of energy, (b) The graph shows how energy 
gap goes to zero. The continuous line describes a second order 
or continuous transition (critical) and the dashed line shows 
the first order transition. The two are distinguished by the 
behaviour of the slope at u — u c . 



probability density does decay to zero but like an unbound 
state it is not normalizable. In higher dimensions, the con- 
dition for a minimal strength of the potential for a bound 
state remains true, but the state corresponding to E = 
becomes normalizable as it should be for a bound state. 
So we see that this bound to unbound transition for a po- 
tential well has different nature in different dimensions. 
In general, (i) for d < 2 there is no such transition as 
E = requires Vq = 0, though there are remnants of the 
transition as Vq — > 0, (ii) for 2 < d < 4, the transition is 
continuous (critical) — the bound state becomes unbound 
through a non-normalizablc critical state as we change u, 
and, (in) for d > 4, the bound state remains normalizable 
up to and including E = 0, and becomes unbound as u is 
decreased further, thus making the transition first-order. 
This depicts a QPT and the case of a potential well gives 
a simple example of a quantum critical point for 2 < d < 4 
with diverging length scales. 

The ground state energy, for u close to u c , is the gap 
A in the spectrum. A quantum phase transition is char- 
acterized by a vanishing gap. A discontinuity of the first 
derivative dA/du signals a first order transition, otherwise 
it is critical or continuous, as shown in Fig. [TJ One may 
define characteristic time and length scales 



£|l = ftA" 1 , and = H/v2mA, 



(2) 



both of which diverge as A — > 0, with £ii ~ z (= 2 in 
this case) being the dynamic exponent. One may compare 
with the classical ground state to see the importance of 
quantum (zero-point) fluctuations and the importance of 
time or dynamics in quantum phase transitions. A path- 
integral interpretation of these scales, useful for the DNA 
mapping, is given below. 

Let us now consider the ground state of two dissimilar 
particles interacting via a central potential V(|rx— r2|) of 
the type of Eq. [TJ with r, denoting the co-ordinate of the 
i-th particle. The existence of diverging length scales and 
scaling behavior around u = u c justifies the dissociation of 
the bound pair to be a QCP or a QPT depending on the 
dimensions they are in. The criticality is described by the 
exponents for the diverging length scales and the energy, 
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Fig. 2: Path integral representation in the x-t plane, (a) A rel- 
ative coordinate path for two particles in one-dimension. The 
solid portions represent the classical bound state, i.e., inside 
the well (B), and the dashed portions represent the unbound 
(U) state in the classically forbidden region, (b) Correspond- 
ing path representation of two quantum particles with time, 
though intersections of paths are not shown explicitly. It is 
also a configuration of two classical Gaussian polymers inter- 
acting at the same contour length as for DNA base pairing, the 
t-axis representing the contour length (z) of the polymers. The 
dotted lines are the melted bubbles whose partition functions 
are characterized by the reunion exponent This description 
holds for any general d. 



and £j 



u — u c 



(3) 



with 



zu± = 



- 1), for 1 < * < 2, (4) 
1, for * > 2, (5) 



which involve (i) z the dynamic exponent, and (ii) & uni- 
versal exponent known as the reunion exponent for 
polymers [91-111] . For the short range interaction problem, 
"J = d/2, as for random walkers, from which the specialty 
of d = 4 is apparent. 

In a quantum bound state a particle can tunnel through 
the potential. In a path integral approach the particle 
does a sizable excursion in the classically forbidden region 
outside the interaction well, sooner or later returning to 
the well (see Fig. [5]). That the two particles will even- 
tually be close-by to form a bound state is the source of 
entanglement while the excursions produce spreads of the 
trajectories in space and time. These spreads give the two 
relevant length scales £|| , The large width of the bound 
state wave function near the QCP ensures the mutual in- 
fluence of the particles even if far away from each other 
(r ^> a) so that the reduced density matrix for one particle 
still carries the signature of the entanglement and the crit- 
icality. For this bipartite system, we are interested in the 
"particle-partitioning entanglement" [13] . This makes the 
von Neumann entropy a valuable quantity for the transi- 
tion which reads, 



S = — Tr p lnp, 



(6) 
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where p is the reduced density matrix for the ground state 



-A„d-A 



p(r ls ri) =Tr 2 g(l,2) = / d d r 2 (r^r^)^ 



(7) 



obtained from the two particle density matrix 
g(l, 2) =| tp)(fp | by integrating out (or tracing out) 
particle 2. In Eq. [SJ we shall introduce some pre-chosen 
length scale to make the argument of log dimcnsionless. 
If, with rrii , denoting the mass and the position of the 
ith particle, the full ground state wave-function (including 
the center of mass (CM)) is 



V>(ri,r 2 ) = $ 



miri + m 2 r 2 
m i + m 2 



<p(ri - r 2 ), 



(8) 



where $ is CM wave function (plane waves) and ip is the 
wave function in relative coordinate (the relative wave- 
function), then 



p(ri,ri)= / d d r 2 ^(r 1; r 2 )^(ri,r 2 ). 



(9) 



Although the center of mass and the relative parts are not 
entangled, the two particles are entangled. The lack of 
knowledge of the state of one particle is the source of a 
nonzero entropy associated with the reduced density ma- 
trix PUS]. 

The translational invariance of the interaction guaran- 
tees that the reduced density matrix p(r,r') = p(r — r') 
has exp(— iq ■ r) as the eigenvector, 



Jd d v'p{r-v')e-^ r ' = p(q)e~ 



with the eigenvalue 



/5(q) = 



Kp 

m 2 



(10) 



(11) 



K being CM wave vector and 0(q) the normalized mo- 
mentum space wave function, the Fourier transform of the 
relative wave-function <p(r) in Eq. [8l Since the entropy 
involves an integral over the whole range of q, it is inde- 
pendent of the CM wave-vector, an expected consequence 
of Galilean invariance. Therefore, without any loss of gen- 
erality, we choose | K |= 0. The eigenvalues constituting 
the "entanglement spectrum" can be written in a scaling 
form 

\<b(q)\ 2 = K- d F(q//t J o«), (12) 

where k 2 = 2p\E\/h 2 = £j 2 , p being the reduced mass. 
Eq. [12] satisfies Tr p = 1. In the critical regime (also called 
the "unitarity limit"), are — > 0, if the scaling function be- 
haves smoothly, then 

F(q,are)^F(q,0) = /(q), (q = a/re) (13) 

which we find to be true for d < 4. For d > 4, we find 
that F(q, an) for ak^O behaves in a singular fashion as 



F(x,y)^y d - 4 /(x) 



(14) 



Here 



so that the prefactor in Eq. [T2] becomes K 
f represents a generic function. By using these limiting 
forms, we find the entanglement entropy to be 



S = P Irian + Co, 
P = min(d, 4) , and cq 



(15a) 

Jd d x f(x)\nf(x). (15b) 



The last statement can be verified by direct computation 
of the momentum distribution function of the relative mo- 
tion in d-dimensions. There are further log-corrections at 
d = 2 and d = 4 which we do not discuss here. 

To motivate Eq. I15al let us consider a few examples. 
Consider the quantum problem of two particles interacting 
via a delta- function potential in one dimension: V(x) = 
—vq8(x). By using the center of mass and the relative 
coordinate wave-function, we write the wave function as 



C e 



iKfil 



-k\xi— x 2 \ 



(16) 



which is translationally invariant. Here K is the CM wave 
vector, k = CJ 1 , and C is the normalization constant. The 
reduced density matrix for particle 1 is then 



p(x, x 1 ) 



C 2 



- {iK [i/rrt2 +k)\x' —x\ 



1 + k\x — X 



having eigenvalues (Eq. [TT]) 
2 1 1 

/5(q) 



ttk(1 + 5 2 ) 2 



(K = 0), 



(17) 



(18) 



which is of the form Eq. [T3]with /(q) ~ (1 + q 2 )~ 2 . By 
introducing an arbitrarily chosen well strength v or a scale 
a = h 2 /2pv in Eq. [SI the entanglement entropy is found 
to be of the form of Eq. I15al with 



P = l, and c = ln8?r - 2. 



(19) 



For k — > 0, p(q) — > S(q) with S = 0. There is a difference 
between n — > and k = 0. 

For a one-dimensional problem with the potential of Eq. 
[TJ one can go over to the delta function potential problem 
by taking a — > keeping Vqcl = vq constant to get the 
same In k behaviour as in Eq. 1191 

We then check for a 3-dimensional potential well, Eq. 
[T] The relative wave-function (I = 0) for this potential is 



ip{r) 




r < a 
r > a, 



(20) 



with k and k as defined earlier and constants A, B deter- 
mined in the usual way of continuity of the wave function 
and its derivative. A direct Fourier transformation of ip(r) 
has been used to numerically compute the entanglement 
entropy. To derive an analytical formula, we note that the 
dominant contribution in /5(q) in the limit an — > comes 
from the outer part. In this approximation we get 



i i 

tj.3 /T7-2 



1 



V(q). 



(21) 
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This p(q) satisfies the normalization condition 
J d 3 q |5(q) = 1. Thus for the 3D potential well in- 
teraction, the entanglement entropy is of the form of Eq. 
II Sal with 



P = 3, and c = 2(1 + log(47r)) 7.06205. 



(22) 



Exact numerical computations of von Neumann entropy 
for d = 3 arc done by using MATHEMATICA. For a given 
k with a = 1, we determine Vb, the depth of the well and then 
the matching conditions and the Fourier transform were used 
to obtain the entanglement spectra. The entanglement entropy 
is then obtained by a numerical integration. The results are 
shown in a log-linear S vs. ft plot in Fig [3] which also shows 
the line obtained from Eq. I15al and Eq |22l It shows that S is 
negative for small k and that it has linear In k dependence. The 
approximations show that the entropy is determined mainly by 
the outer part of the wave-function. 




Fig. 3: Plot of S vs. In k with a = 1. The circles are the 
numerical values and the straight line is the predicted line S = 
3 In k + 7.06, Eq. [H 

To generalize the result for any dimension we carried out the 
calculation for general d. The density matrix, solely from the 
outer part, is expected to be of the form /(q) ~ (1 + q 2 )~ 2 
as in previous cases but then there is a divergence problem for 
normalization for d > 4. Since we want Tr p = 1, an ultraviolet 
cutoff is required. This makes an an important variable even 
in the limit an — >• 0. The specialty of d — 4 is now evident. 

The radial wave function R(r) (1 = state as the ground 
state) is, 



R(r) 



where r 



A W 2 J| £/2 |(fcr) 
B r^H^AiKr) 



for r < a 
for r > a. 



(23) 



2 — d. A, B determine the normalization and 
matching of the inner and the outer solutions. Here J and 
ffW are the Bessel and the Hankel function of the first 
kind. The continuity of the wave function at r 



AJ\2-d\{k r ) 



B(na) 



a gives 
(24) 



under the condition k — > and ka — > k c = tt/2. Eq|3] 
follows from Eql24l the matching of log derivative and 
the Bessel function identities. By using the normalization 
condition and Eq. [24l we get 



B 



,\2-d\/2 



for d < 4 
for d > 4. 



In the same k — > limit, with outer part dominance, 



4>{q = 0) « Bk 



(25) 



which gives 



p(q) = |0(q)| 2 = B 2 k 



-(2+d) 



/(q) 



for d < 4 



for d > 4. 



(26) 



So the von Neumann entropy is of the form Eq. (|15a[) with 
P = 4 for d > 4. 

In terms of the deviation from the critical point, the 
entropy is 



In lu- 



for * < 2. 



(27) 



For the case in hand, \P = d/2. The form of Eq. |2"T1 brings 
out the universal behavior of the entropy and has validity 
for potentials different from Eq. [1] like e.g. scale- free 1/r 2 
potential JT2] . All the details of the interaction go in the 
universal exponents z and \& '. The entropy diverges at the 
critical point and, is negative. 

DNA connection. — We show the connection of the 
quantum entanglement entropy to the entropy of bubbles 
in DNA melting. Under an imaginary time transforma- 
tion, the path integral formulation of the quantum prob- 
lem is analogous to a classical statistical mechanical sys- 
tem of polymers used in the context of melting of DNA 

latEHDHi. 

Let us consider a DNA whose two strands are two 
Gaussian polymers in <i-dimcnsions and index the points 
(monomers) by the contour length z measured from one 
end. The native base pairing of a DNA requires that a 
monomer at index z on one strand interacts with a point 
on the other strand with the same index z. This is the 
Poland-Scheraga type model [TT] for DNA melting. By 
using one extra coordinate for the sequence or the length 
of the polymers, we get directed polymers in d + 1 dimen- 
sions like paths in path integrals, as shown in Fig [5J In 
this representation the base pairing interaction maps onto 
the same time interaction of the quantum system, time 
playing the role of the base pair index. The DNA parti- 
tion function as a sum over all polymer configurations is 
equivalent to the sum over all paths in quantum mechan- 
ics. The DNA Boltzmann factor cxp(— j3H) with j3 as the 
inverse temperature and H the Hamiltonian for two chains 
of elastic constants Kj as 



/3H-- 



N 



3=1,2 



dz 



-V(v 1 (z)~r 2 (z)) 



dz, (28) 



corresponds to the factor exp(iS/H) for path integrals with 
S the classical action of two interacting particles under 
z — > it. This makes the Green function or the propagator 
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G(x\, x^i t\xi, x' 2 , 0) equivalent to the partition function 
Z(xi,X2,N\x' 1 ,x' 2 ,0), (N — > ir). Here Xj,x'j are the co- 
ordinates of the j-th strand end-points at and at length 
N. The free energy per unit length of DNA for N — s- oo 
is the ground state energy of the quantum problem. 

The short range base-pairing potential can be taken to 
be a contact potential or a well of Eq. (JTJ . Then the pic- 
ture of return of the quantum particles within the range 
of interaction after excursions outside the well gives the 
equivalent picture of polymers with broken base pairs hav- 
ing excursion away from binding and eventually coming 
back to the well to form pairs. This excursion swells the 
polymer and creates bubbles along the length of the DNA. 
Thermal energy opens up bubbles in the bound state of 
DNA. The entropy of a bubble of length N is determined 
by the reunion partition function of two polymers starting 
together and reuniting again at N, which for large N, has 
the form Q(N)=N~^e Na °, or the entropy 

S = lnfi(iV) = Na -*lnJV, (29) 

in units of the Boltzmann constant feg = 1. Eq. 1291 shows 
that do is the bubble entropy per unit length that survives 
in the thermodynamic limit. However, the power law N- 
dependence which gives the negative sub-extensive part 
of the entropy is essential for the transition and also for 
the bound state. The reunion exponent ^ determines the 
universality class of the binding-unbinding transition and 
there is a melting transition if and only if 'f? > 1 . See Ref. 
[H] for details. 

The one-dimensionality of the chains requires an alter- 
nating arrangement of bound regions and bubbles as in 
Fig. [5J The arrangement allows one to write the parti- 
tion function, after Laplace transform with respect to the 
length (i.e. in the grand canonical ensemble) [9], as 

G(x,y;a) = G Q {x; s)G(0, s)G Q {y; s) 
G Q {x; s)G Q (y; s)G B (s,u) 
l-G u (.s,a )G B ( S , U ) ■ {6U) 

Here x = {xi,X2},y = {x^x^}, G is the Laplace trans- 
formed partition function of two polymers tied at one end 
and open at the other, called the survival partition func- 
tion, and G(0, s) is the total partition function with two 
ends bound. In G Q , the tied point is to be integrated over 
keeping the set x or y fixed. G(0, s) can be written as 
a sum of a geometrical series (see Fig. ^ involving the 
partition functions of the bound parts and the bubbles, 
G B (s,w) and G u (s, <jq). The free energy comes from the 
singularity of G(x,y,s) which is either s = ao = or at 
s = so for which 

G JJ (s,* )G B (s,u) = l, (31) 

with (To — 0, so satisfies Eq. ((3]). 

Near the nontrivial singularity, a pole at s = so, the 
form of G(x, y, s) resembles the Green function in the en- 
ergy eigenfunction expansion as 

(y\tp)(ijj\x) 



with ground state dominance. From the equivalence be- 
tween DNA model and the quantum problem, we identify 
the density matrix as 

p(x,y) ~ G (x;so)G (y;s )/G v (s ), 

so that the entanglement entropy would behave like 
S ~ In G u (soi Co). By using Gaussian distributions for 
Gaussian polymers (i.e. random walkers), one recovers 
Eq. HU 

To get the behaviour of lnG u , we employ a finite- 
size scaling analysis. The phase transition in the poly- 
meric system occurs in the N — > oo limit so that a fi- 
nite N acts as a finite size scale both for DNA and in 
the quantum problem. The finite size scaling variable is 
N/£j_ so that the entanglement entropy is proportional 
to -z\n£j_ - ^ In | u - u c | (see Eq. [27J). The dif- 
ference in the amplitude occurs because of the different 
normalization used for polymers and quantum problems. 
The point to note is that the entanglement entropy in the 
quantum problem comes from the universal non-extensive 
part of the entropy of the bubbles. Since the full entan- 
glement spectrum is known, it is also possible to compute 
the Rcnyi entropy [6] . We recover in the appropriate limit 
the result quoted in Eq. I15al In the DNA interpretation, 
the Renyi entropy would come from many circular single 
strands (replicas) pairing with a large single strand, re- 
sembling the rolling circle replication of viruses. Details 
will be discussed elsewhere. 

Discussion. — A negative entropy is counter-intuitive 
when one has the third law of thermodynamics in the back 
of one's mind, though exceptions are known; e.g. neg- 
ative entropy is found for perfect gases at low tempera- 
tures or as a corollary of the classical equipartition theo- 
rem. One can see the same feature by writing the reduced 
density matrix in terms of an entanglement Hamiltonian, 
p oc cxp(— /3-ffent)) m a form reminiscent of a Boltzmann 
factor. The diagonal form in Eq. [2~T1 shows 

f3H cnt = 2ln(l + q 2 /K 2 )^2q 2 /K 2 , (for small q), (33) 

which is like a classical d-dimensional oscillator in g-space, 
with k 2 as the effective temperature. A direct calcula- 
tion or use of the classical equipartition theorem now tells 
us that the entropy has d In k behaviour as in Eq. I15al 
We believe this to be a generic feature whenever the en- 
tanglement Hamiltonian is gapless. Another way to see 
this emergence of In k in entropy is to compare with the 
DNA problem. The equivalent classical DNA model also 
has a negative diverging part of entropy but that sub- 
extensive part vanishes in the thermodynamic limit of the 
entropy per unit length. In the quantum case, the equiv- 
alent limit has no such advantage in finding the entropy 
because demanding extensivity in time direction is mean- 
ingless. Hence the negatively diverging term is inevitable 
near criticality. 

In this paper we show that the quantum entanglement 
entropy near the bound-unbound transition of two inter- 



p-5 



Poulomi Sadhukhan and Somendra M. Bhattacharjee 



acting particles comes out to be negative, and it diverges [15] Rio L., Aberg J., Renner R., Dahlsten O. and Ve- 
at the QCP. Using the equivalent classical statistical me- dral V., Nature, 474 (2011) 61. 

chanical system of DNA near the melting transition we 
show that the negativity of the entanglement entropy is a 
necessity and is essential for the phase transition. The co- 
efficient of the logarithmic term contains the information 
of the interaction and the universal behaviour of the phase 
transition. The coefficient is shown to be related to the re- 
union exponent of vicious walkers. This is the first time in 
the context of quantum entanglement that the negative en- 
tropy is found by explicit calculation. We argue that this 
log divergence in the quantum case and the sub-extensive 
part in the DNA problem are linked by finite size scaling 
near the critical point. From the renormalization group 
(RG) approach for the DNA melting problem [TU1I12] . one 
may infer that the entanglement entropy increases along 
the RG flow, since the critical point corresponds to the 
unstable fixed point. It has been argued recently that en- 
tanglement can be used to produce negative entropy |15j . 
The information theoretical meaning of the negative en- 
tropy in our case is not very clear. Our speculation is 
that the negative entropy is the norm, not an exception 
near a quantum binding-unbinding transition. We feel sig- 
natures of negative entropy might be detectable in cold 
atoms where interactions can be tuned to the unitarity 
limit. If one can harness the negative entropy, one may 
cool a system or a computer and possibly may overcome 
the obstacle to circuit miniaturization. 
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